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Abstract 

This paper deals with monic orthogonal polynomials generated by 
a Geronimns canonical spectral transformation of the Laguerre classical 
measure, l.e., 

—1— x°‘e~^dx + NS{x — c), 

X — c 

for X € [0,oo), a > —1, a free parameter N £ R+ and a shift c < 0. 
We analyze the asymptotic behavior (both strong and relative to classi¬ 
cal Laguerre polynomials) of these orthogonal polynomials as n tends to 
infinity. 

AMS Subject Classification: 33C45, 41A60, 33C15 
Keywords: Orthogonal polynomials, Canonical spectral transforma¬ 
tions of measures. Asymptotic analysis. Hypergeometric functions. 


1 Introduction 

Let /r be a positive Borel measure supported on a finite or infinite interval 
E = supp(/i), such that the convex hull verifies Cq{E) = [a, b] C M. In the last 
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years several authors have studied the so called canonical spectral transforma¬ 
tions of /i, which are a way to construct new families of orthogonal polynomials 
from a perturbed version of p. They have been studied from several points of 
view, including the corresponding Jacobi matrices (see 0, m) or the Stieltjes 
functions associated with such a kind of transformations (see [5H] among others). 

Let us introduce the sequences of monic orthogonal polynomials (SMOP in 
the sequel) associated with one of the aforesaid canonical transformations, called 
the Geronimus canonical transformation on the real line. The basic Geronimus 
perturbation of /i is defined as 

- dp,{x) + NS(x — c), 

X — c 

where N G M+, 5{x — c) is the Dirac delta function located at a; = c, and the 
shift of the perturbation verifies c ^ E. Observe that it is given simultaneously 
by a rational modification of /i by a positive linear polynomial whose real zero 
c is the point of transformation (also known as the shift of the transformation) 
jointly with the addition of a Dirac mass at the point of transformation as well. 

This transformation was introduced by Geronimus in his pioneer works m 
and m on procedures of constructing new sequences of orthogonal polynomi¬ 
als from other known sequences families, and it was also studied by Shohat 
in a different scheme involving mechanical quadratures (see [23]). Years later, 
Maroni (see HHI) returned to the problem and gave a first expression of the 
Geronimus perturbed orthogonal polynomials in terms of so called co-recursive 
polynomials of the classical orthogonal polynomials. More recently, Bueno and 
Marcellan reinterpreted this perturbation in the framework of the discrete Dar- 
boux transformations (see |4|). In |3] the authors present a new computational 
algorithm for computing the Geronimus transformation with large shifts. In 
|2j the authors provide sharp limits (and the speed of convergence to them) of 
the zeros of the Geronimus perturbed SMOP, and also, when p is semi-classical 
they obtain the corresponding electrostatic model for the zeros of the Geronimus 
perturbed SMOP, showing that they are the electrostatic equilibrium points of 
positive unit charges interacting according to a logarithmic potential under the 
action of an external field. In m the authors extend the standard Geronimus 
transformation to a cubic case. [10] provides a new revision of the Geronimus 
transformation in terms of symmetric bilinear forms in order to include certain 
Sobolev and Sobolev-type orthogonal polynomials into the scheme of Darboux 
transformations. Finally, |9| deals with multiple Geronimus transformations 
and show that they lead to discrete (non-diagonal) Sobolev type inner prod¬ 
ucts, and it is shown that every discrete Sobolev inner product can be obtained 
as a multiple Geronimus transformation. 

In view of the foregoing, this transformation has been extensively studied 
in the literature, mainly in analytic and algebraic frameworks. However, to the 
best of our knowledge, the asymptotic properties of the family of orthogonal 
polynomials as u —>■ oo have not been studied in detail, save for the particular 
case when N = 0 and the perturbed measure is the Laguerre classical measure 
(see in]). 
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2 Laguerre polynomials and functions of the sec¬ 
ond kind 


The classical Laguerre polynomials L“(x) are defined as the polynomials or¬ 
thogonal with respect to the inner product 

poo 

{p,q)a= / p{x)qix)x°^e~^dx, a >-1, p,q€¥, 

Jo 

see, among others [S] or |22j . 

In order to hx notation, we denote by L^(x) the monic Laguerre polynomials, 
so L“(x) = x” -I- .... These monic polynomials are connected to standard 
Laguerre polynomials L“(x) by the formula 

L“(x) = (-l)"n!L“(x), n > 0. (1) 


They satisfy a three term recurrence relation that we write in the following form 
xL“(x) = L“+i(x) -b/3„L“(x) -b 7 „L“_i(x), (2) 


where 


= 2n-I-Of-I-1, jn = n{n + a), 


(3) 


and we have initial data Lq (x) = 1 and Lf{x) = x — a — 1. We will also make 
use of the L^([0,oo)) norm of the monic Laguerre polynomials. Since 


a||2 _ r(n -I- Of -I- 1) 

r(n-bl) 


we have 

||L“||2 =r(n + a + l)r(n + l). 

A second (independent) solution of the recurrence relation (U) is the function 
of the second kind, obtained via a Stieltjes transform of the Laguerre polyno¬ 
mials ^ 

Kiz) = / (4) 

Jo z 

which is an analytic function for z G C\ [0, oo). Using the Rodrigues formula for 
Laguerre polynomials [U §18.5(ii)] and a standard integral representation, see 
[U Eq. 13.4.4] for instance, it is possible to write F^{z) in terms of the confluent 
hypergeometric function of the second kind, or Kummer U function: 

E“(z) = (-l)"n!r(n + a + l)U{n + 1,1 - a, (5) 


with plus sign if —tt < argz < 0 and minus sign if 0 < argz < tt. This 
representation will be a key element for all the asymptotic analysis later on 
in this manuscript. For more information about the confluent hypergeometric 
functions, we refer the reader for instance to [1] Chapter 13]. 
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Let us introduce the following inner product in the linear space P of poly¬ 
nomials with real coefficients 


pOO 

{p,q)vN = / P{x)q{x) 


X — c 


-dx + N6{x — c), 


( 6 ) 


where a > —1, and N > Q, and c S (— 00 , 0). Namely, we deal with a measure 
that consists of an absolutely continuous part, which is a rational perturbation 
of the Laguerre weight on [0, -l-c»), plus a Dirac delta located at point x = c: 


dvN{x) 


X — c 


dx + N6{x 


c). 


Equivalently, we will say that is a Geronimus perturbation of the standard 
Laguerre measure (see mm and the references therein), and we will denote by 
{x) the monic orthogonal polynomials with respect to (|6]). 

As explained in [5] (and the references therein), the Laguerre-Geronimus 
monic orthogonal polynomials can be written in terms of the monic Laguerre 
OPs using the following simple connection formula: 

Qr’^ix) = LZix)+A^LZ_dx), (7) 

where the coefficient depends on n, a, c and N. It can be written as 


A 


N 

n 


r(n -I- Q!)r(n) 

L“_i(c)F“_i(c)-iVL“_i(c)2 


( 8 ) 


where we have defined, for c G (— 00 , 0), 


7r„(c) 



(9) 


We observe that in Proposition 1 and Remark 1 in |2], different expressions 
for A^ are given, but using identities for the Laguerre polynomials it is not 
difficult to obtain ([5|) from those. 

A particular case is given by A^ = 0, where the Dirac delta is not present in 
the perturbed measure. Then A^ = —r„_i(c), where 


^n(c) = , ce(-OO,0), (10) 

with F^^{c) = 1, analogously to |12l §2.4.4]. 

The aim of this paper is to obtain strong and relative asymptotics of the 
sequence of OPs Q“’°’^(.z), for large n and other parameters fixed. The sim¬ 
plicity of the connection formula ([7]) makes it a very attractive identity to use 
in conjunction with classical asymptotic approximations for Laguerre polyno¬ 
mials (such as Perron, Fejer or Mehler-Heine expansions, see [2^ §8.22]), in 
order to obtain the corresponding result for the Laguerre-Geronimus OPs. The 
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only element that is missing so far in the literature is a study of the asymptotic 
behavior of the coefficient . We observe that because of ([8]), depends on 
the ratios 7r„(c) and rn(c). 

The structure of the paper is as follows: in Section [3] we obtain large n 
asymptotic expansions for 7r„_i(z) and rn-i(z) in the complex plane, which 
lead to asymptotic approximations for A^ in Section 01 Putting together this 
result and the connection formula we obtain in Section the strong and 
relative asymptotics for {z) in the complex plane. Finally, Section [3] is 

devoted to the study of the coefficients in the three term recurrence relation 
satisfied by and Section [7] gives a hypergeometric representation for 

the OPs that could be of independent interest. 


3 Asymptotic expansions for TCn-i and r„_i 


The ratios iTn-iiz) and rn-i{z) could in principle be studied using standard 
techniques for the asymptotic behavior of solutions of three-term recurrence 
relations, such as the Perron theorem, see for instance |151 §4.3]. However, since 
the recurrence coefficients in m satisfy /3„ ^ 2n and 7 „ ^ as n —>■ oo, the 
theorem is inconclusive about the existence of minimal and dominant solutions, 
and it does not give detailed asymptotic information about the behavior of ratios 
of solutions. We refer the reader to [T] Section 4] for more details. 

In this paper we work with strong asymptotics of the Laguerre polynomials 
and functions of the second kind directly. For 2 away from [0, oo), the strong 
asymptotics for the L‘^{z) can be obtained from the classical expansion due to 
Perron, see for instance [121 Theorem 8.22.3]: 

L“(z) = (l + , (11) 

which is valid for fixed a > —1 and z G C \ [0,oo). The fractional powers are 
assumed to take their principal values, with phase between —tt and tt. In (51 
Theorem 3] higher terms in this asymptotic expansion have been obtained, using 
a related expansion for confluent hypergeometric functions due to Buchholz, see 
also [16| and references therein. The ratio asymptotics is given in |6] as well: 




1 + 
1 + 



2a — 2^; — 1 
4(n — 1) 





2a — 22 — 1 
4n 


-kO(n-3/2), 


as n —>■ 00 , for fixed a and z G C \ [0,oo). Therefore, as a direct consequence, 
for the monic polynomials we have 


1(2) = 






= —n — 


2z — 2a + 1 


(12) 
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as n —>■ oo. 

Regarding the asymptotic behavior of the functions of the second kind, we 
present the following result: 


Proposition 1 Given fixed c € (—oo,0) and a > —1, the functions of the 
second kind F°(z), defined by n, satisfy 


F^(z) = (-l)-^(-z)f-3 e-"/2-2^^r(n - 


l)r 


eo 




(13) 


The expansion is valid for bounded z S C \ [0, oo), with principal values of the 
power functions. The first few coefficients Cj are 


eo{a,z) = 1, 

, , 12a2 _ 3 _ 24z(1 - a )- 

ei(a,z) =--, 

, , 16z4 + 192(1 - a)z3 + 24(200^ - 48a + 13)z2 

=-4608- 

144(a - l)(2a + l)(2a + 3)z + 9(40^ - l)(4a2 - 9) 

4608 ■ 

We observe that the leading term in this expansion is consistent with the 
results in mms], bearing in mind that we are working with monic polynomials. 
We also note that the exponential factor is erroneously corrected in Proposition 
3.2 (a) in the first reference. 

As a direct consequence, and using symbolic computation, we have an asymp¬ 
totic expansion for the ratio of consecutive functions of the second kind: 


Proposition 2 As n —> oo, the ratio asymptotics of the Laguerre functions of 
the second kind is given by 


rn-i{z) 


Fuiz) 


, - 2z 2a “t” 1 , _ 1 jr. 

-n + ^-zn H- - - l-0(n ' ), 


where a > —1, c € (—oo, 0) and z G C \ [0, oo). 


We remark that division by Fff_i{z) is allowed for z away from the positive 
real axis, bearing in mind ® and the fact that the Kummer function U (a, &, z) 
does not have zeros for | arg z| < tt if a is positive, see [TJ §13.9]. 

In order to get the previous results, there are at least two possibilities: use 
the expression of F^{c) in terms of Kummer functions, see or deduce its 
asymptotic behavior from the Deift-Zhou steepest descent method applied to 
the corresponding Riemann-Hilbert problem, see the work of Vanlessen [26| and 
the monograph by Deift [H] . In the sequel, we elaborate on the first approach. 
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Following the ideas exposed in [21], we use the integral representation for 
the Kummer [/-function: 

1 

U{a, b,z) = ^ 

r(a) Jo 

which holds for Re a, Re z > 0, see also [TJ 13.4.4]. The transformation //(l-|-t) = 
e~'^ gives 

r°° , 

Uia, 6, z) = ._ / e-“"-^/V-V(r)dr, 

r(a) Jo 

where 

/(-) =(tz^ 

The function /(r) is analytic for |r| < 2tt, and therefore it admits a power 
series expansion around the origin of the form 

OO 

/('t) = X! Jyn{h,z)T^. (14) 

m—0 

Integration term by term, invoking the classical Watson lemma, [201 [25] , 
gives an expansion of the [/-function of the form 

M-l 

U{a,b,z)= ^ dm{b,z)(j)m{a,b,z) + RMia,b,z), 

m—0 

where the asymptotic sequence is 

(j)mia,b,z) = ^:^—(-) ' R:^+i_h(2(az)^/^), 

i (a) \a/ 

in terms of modified Bessel functions, using the fact that 

R,(2(zC)1/2) = ^ (Q ' e— 

valid for Rez,ReC > 0. Also, we have RM{a,b,z) = 0{(j)M{a,b, z)) as a —>■ oo, 
uniformly with respect to z in compact sets in z > 0 and uniformly with respect 
to b in compact sets of R. We assume, following |24l §2.1] that M is large 
enough, in particular M > b. 

If we replace b = 1 — a, symbolic computation gives 



do{a,z) = 1, 


di{a,z) 

d2{a,z) 

d3{a,z) 


6(1 — a) — z 
12 ’ 

z2 - 12(1 - a)z + 12{a - l)(3a - 2) 

288 ’ 

—5z^ — 90(1 — a)z^ — 36(15a^ -I- 25q; -I- 8)z — 1080Q;(a 

51840 


1)2 
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and so on, for the coefficients in the series expansion m- Replacing a by n + 1 
in the asymptotic sequence, we obtain 


0m(n,a,z) = R:„+„(2((n+(15) 

n\ \n + I / 

As a consequence, using ([S]) and (ITCl) . we have 

F“(z) = (-l)”n! r(n + a + 1) C/(n + 1,1 - a, 

= 2e“^/^(-l)"r(n + a + 1) [SAiin, a, + Rnin, a, ze=*"”)] , 

as n —>■ oo, where 


M-l 

SM{n,a,ze^^") = ^ dm{a,-z) 
m—0 



m-\-a 

if^+„(2((n + l)ze±-*)i/"). 


and a, is the remainder. Once again, we take plus sign if —tt < 

arg z < 0 and minus sign if 0 < arg z < n. 

It is possible to re-expand this asymptotic series in inverse powers of n: using 
the asymptotics of the modified Bessel functions for large values of the argument 
and fixed order z/, see [3 10.40.2], we have 


K, 


e=o 


ai{v) 


00 , 


where ao(z^) = 1 and for £ > 1, 

- 1 )( 4 z .2 _ 32)... _ (2£ _ 1)2) 


ai{v) = 




So, if we denote s„ = 2((n + l)ze^’^®)^/^, we have 


K, 


m+a 


(Sn) 



ae{m + a) 


Note that s„ = 2'/nze^^ {1 + 0{n~^)) as n —>■ 00 . Assembling all the previous 
results and expanding in inverse powers of n, we arrive at Proposition [TJ Propo¬ 
sition [2] follows then from this result, using symbolic computation to manipulate 
the asymptotic expansions. 


4 Asymptotic behavior of 

The main result of this section is the following: 









Proposition 3 Let c G (— 00 , 0), a > —1 and N > 0 be fixed parameters, then 
= n + \/—cn + —— ^ — - + 0(n~^^^), n —>■ 00 , (16) 

if N > 0, and 


K. = =n - 


2a — 2c — 1 


+ 0(n-i/2), n 


(17) 


Proof. We deduce this result from the asymptotic expansions derived before: 
combining dm with m, we obtain 


L“_i(c)F“_i(c) = 


1 


2^=7 


:r(n)r(n +a) (1 + 0(71-1)). 


as 71 —>■ 00 . It can be checked that the term of order 0(n 1 /^), which is to be 
expected here, is actually equal to 0. Also, 

L“_i(c) 2 = (l + 0{n-^/^)) , 


2^/^cn 


: r(77)r(77 + a)el^^ (1 + 0(77-1/2)^ , 


2^c 

using the fact that r(n + Q:)/r(77) = 77 “(1 + 0{n~^y), as 77 —>■ 00 , and the 
notation 

e“(-c)-“ 


Lt(y^C - 


27r 


Putting everything together and using the results in [^, we have 

2y/—cn 


A^ = - 


Lu{c) 


1 + 0(77-1) - AfO„,,e4V=^(l + 0(77-1/2)) A“_i(c) ’ 

2^/—cn 

1 + 0(77-1) - AfO„.ee4^/^(l + 0 ( 77 - 1 /^ 

2a — 2c — 1 


+ 77 + y/—cn + 


+ 0(77-1/2), 


It is important to observe that if > 0, the first term in the previous 
sum is exponentially small in 77 , so it does not contribute to the final result. 
However, if A^ = 0 the first term does contribute, since the exponential term 
is not present, and then we have the difference in sign in the subleading term 
given in Proposition [3l ■ 


5 Asymptotics for 

Using the estimates for (HD) and (El), we describe first the strong asymp¬ 
totics for for 0 away from the interval of orthogonality (outer asymp¬ 

totics): 
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Proposition 4 Given fixed values of N > 0, a > —1, c € (—oo,0) and z G 
C\ [0, c»), as n ^ oo, the monic polynomials (z) verify the following 

strong asymptotics 




{V^ztV- 

(18) 

and the following ratio asymptotics with respect to monic Laguerre polynomials: 


^ \/^T \/^ _ (\/^T 
L^iz) Vn 2n 


(19) 


In both cases, the upper sign corresponds to the case iV > 0 and the lower sign 
to N = 0. 


Proof. We rewrite the connection formula ([71) as follows: 


^a,c,N 


(z) 



K-iiz) \ 

K{z) ) 


Kiz), 


( 20 ) 


and we use the information obtained so far. Observe that division by Lf^iz) 
does not cause any problem, since the zeros of this polynomial are contained in 
[0, oo). From (flTl) . we deduce 


K-iiz) 

L^iz) 



I — 2z — 2a 
An 


+ C>(n-3/2) 


n —>■ 00 , 


and combining this with the asymptotic behavior for given by Proposition [3] 
and m, we obtain (fT51) . Then formula (ITIJl) for the relative asymptotic behavior 
is a direct consequence of (1701) . ■ 

We remark that this is consistent with the result in [111 Proposition 3.4 a)], 
taking = 0 and M = 1 . 

Applying the connection formula again, but with the inner asymptotic ex¬ 
pansion for Laguerre polynomials, we can obtain strong asymptotics of [x) 

for X G ( 0 , oo). 

Proposition 5 Given fixed values of N > 0, a > —1, c G (—oo,0) and x in 
compact intervals o/( 0 ,oo), as n ^ oo, the monic polynomials Q“’°’'^(a;) verify 
the following strong asymptotics 


Q 




al2-3l4 x/2 

0 fa.a/ 2 +i /4 sin 6 '“(a;) ± cos 6 »“(a;) -b C>(n ' ) 

( 21 ) 


where the phase function is 


Olf{x) = 2^/nx — 



( 22 ) 


and again the upper sign corresponds to the case N > 0 and the lower sign to 
N = 0. 
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Proof. We rewrite ([7|) as follows: 


QT’^ix) = (23) 

where we have used the following identity for standard Laguerre polynomials: 

Ll_,{z)=LZ{z)-Ll-\z), (24) 


see for example PQ 18.9.13]. Then, we use the classical Fejer formula for Laguerre 
polynomials, see Theorems 8.22.1, 8.22.2], adapted to the monic case: 


Ll{x) = (-l)"n! 


^a/2-l/4ga:/2 


(25) 


valid for x in compact intervals of (0, oo), with phase function (1221) . 
From the asymptotic expansion of h.^ we deduce that 


KN \N I _ 

1-— = T\/—+0(n-^), ^=l±J— + 0{n-^), 

n \ n n \ n 

where the upper sign corresponds to the case N > 0 and the lower sign to 
N = 0. Using this information and (|25L we expand as n —>■ oo in (l23l) . bearing 
in mind that from the definition cos0““^(a:) = — sind“(a:), and we arrive at the 
result. ■ 

Another useful asymptotic behavior that one can find in the literature is 
the Mehler-Heine type formulas. Is very well known that for j S N U {0}, the 
standard Laguerre polynomials verify 

lim + ^ ^ (26) 

n—>cxD 77,“ ^ ' 

uniformly for z in compact subsets of C, see [521 Theorem 8.1.3], where Jq is the 
Bessel function of the first kind, and the square root takes its principal value. 
Using this result, we can prove the following: 


Proposition 6 For given values of N > 0, a > —1, c G (—oo,0) and z in 
compact subsets ofC, the polynomials verify 


lim 

n—¥oc 


(_l)n 


nl 


w 


t-l/2 


— T J a (2 ) 5 


(27) 


where the upper sign corresponds to the case N > 0 and the lower sign to N = 0. 


Proof. In order to prove (|57l) . we start with the connection formula ([7]) for 
monic polynomials. In terms of standard Laguerre polynomials, recall o , we 
have 


(-1)" 

nl 


QZ{z/n) 


Ll{zln)-^Ll 

n 


{zfn) 


= 1 - 


A5 


A^ 

L“(z/n) + ^L: 

n 


a-l 


{z/n), 
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where we have used (IMl) again. Consequently, 

(-1)" Ql{.zln) ^ A _ ^ 
nl n°‘ \ n 

Next, using the asymptotic expansion for A^, we deduce that as n —> oo, 

1 - —=Ta/—+ 0 (n-i), ^=0(n-i), 

n \ n 

where the upper sign corresponds to the case N > 0 and the lower sign to 
N = 0. Then multiplication by in (1^51) and the use of the Mehler-Heine 
asymptotics for Laguerre polynomials (PSI) give the result. ■ 

Again, we note that this is consistent with m Proposition 3.4 b)], taking 
N = 0 and M = 1. 


Lni^ln) , A" L“ ^z/n) 






(28) 


6 Asymptotics for recurrence coefficients 


The moiiic Laguerre polynomials satisfy the three-term recurrence relation (P|). 
and from that we can obtain an analogous recurrence for the perturbed polyno¬ 
mials 

Theorem 1 The polynomials {z) satisfy a three term recurrence relation 

= (^ - Pn)QT'^{z) - lnQTf{z), (29) 


where the coefficients are given by 

= Pn+ ~ ^n+lt 


7n = 


and j5n and 7 „ are given by m- 


A" 

\N 

^^n-1 


ln-1, 


Proof. Let us denote by [3 and 7 the coefficients of the three term recurrence 
relation for the polynomials (^“’'^’^(z): 

xQT^^iz) = (z) + AQ“’^’^(z) + %QTf[z) 

Using the connection formula ([7]) on both sides of the previous equation we get 


xLf,{z) + A^xLf,_,{z) = L“+i(z) + (A(:+i + ^„)L“(z) 

+ {Pn^n + 7n)^n-l(^) + 7n A^L iL “_2 (z) • 

We use the three-term recurrence relation for the Laguerre polynomials Q on 
the left hand side of the previous equation and we obtain 


L“+i(z) + (/3„ + A(:)L“(z) + (7n + AC^/3 „_i)L“_i(z) + 

= L“+i(z) -I- (AO'^;^ -|-/3 „)L“(z) -I- {jdn^n +7n)^“-l(z) + ln^n -1 


n-2(^) 

i“-2(^)- 
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Since the Laguerre polynomials are a basis for the space of polynomials, we 
obtain the following equations 


Pn + + /3„, (30) 

7n + Pn-1 = PnA^ + 7„, (31) 

A0'7„_i =7„A0Li. (32) 

The theorem follows directly from formulas (1501) and (1321) . ■ 

Remark 2 From substituting the expressions for 7 „ and /3„, we obtain 

the following non-linear recursion for A^: 


aAT _ \N ^ _lrL ^ htzL + 2 

■''■n+1 A TV ^ \N ' 


It follows directly from p3p that 

n+1 


A7i - a? = ^(Af - Af_,) = + 2(n - 1 ) + ir. 


i=3 

This gives the following recursion for An: 

n(n + a) 


An 


Ai 


= -1 

Aiy 


+ 2(n — 1) + + A^. 


(33) 


(34) 


If we define A^ = Qn /Qn-ii then HSfl) becomes 


Qn+i - (^2(n - 1) + ^ + Af ^ pO" + n(n + a)p+i = 0. (35) 

From this result and the asymptotic expansion for A^ , see (1161) and (fTTl) , it 
is straightforward to deduce the large n asymptotic behavior of the recurrence 
coefficients fin and +, in terms of the recurrence coefficients for monic Laguerre 
polynomials, fin and 7 „, given by 

Corollary 1 As n ^ oo, for fixed c G (—oo,0) and N > 0, the coefficients 
Pn and 7 „ of the three term recurrence relation for monic Laguerre-Geronimus 
orthogonal polynomials satisfy 

= (l - ^ ^ + 0{n-^)^ +, 

+ = (l + - T + 0{n-^)^ 7„_i, 

where the upper sign corresponds to the case N > 0 and the lower sign to the 
case N = 0. 
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7 Hypergeometric representation of 

In this section we will derive a representation of the Geronimus perturbed family 
of orthogonal polynomials as hypergeometric functions. For this we need the 
connection formula 0 together with the hypergeometric representation of the 
monic Laguerre polynomials, that can be obtained from [U 18.5.12]: 



(36) 


Theorem 3 The monic polynomials have the following hypergeomet¬ 

ric representation 



where 



Proof. It follows from the connection formula © and (1551) that 




By a straightforward calculation, can be written as 



(37) 


Next, we rewrite the expression in square brackets as 



where 


e: 


•n 


N 


n{n + a — 
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By replacing (1551) into (1571) . we obtain 


(^a,c,N( \ _ (-1)" (« + l)ra {-n)k 

n(a + n) (a + 1) 

k—O 

= {-ir{a + l)r.2Fj 

\ n + a J \ 

This completes the proof of the theorem. ■ 


(l + e^)fe 
, ie^)k 
'-n, l + e^ 
. a + 1, 



Remark 4 The hypergeometric functions 2 F 2 are solutions to a third-order dif¬ 
ferential equation JT\ 16.8.3j. Therefore, Theorem\^ implies that the perturbed 
polynomials are solutions to 

z^y”'-x{x-e!^-a-2)y"-{{e^-n + 2)z-{a + l)e’^)y' + n{e^ + l)y = 0. (39) 

This differential equation can be easily obtained from the holonomic equation for 
the polynomials (see |11 Section f.l]): 

y" + R{z)y' + S{z)y = Q, (40) 


where 


R{z) 

S{z) 


Aff a + 1 

zA((+ {n-A(f){n + a-A(()^ z 
zA(t + {n- A((){n + a) n-l 

z{zA(l + {n - A(l){n + a - Ad)) z 


Indeed, (|39l) is obtained by multiplying the derivative of (go]) by z'^ and adding 
it to (1401) multiplied by —z{z — ed — a — 2) — z'^R{z). 
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